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Computational  Aeroacousties:  An  Overview 


Christopher  K.W.  Tam 

Department  of  Mathematics,  Florida  State  University 
Tallahassee,  FL  32306-45 1 0,  USA 


Abstract 

An  overview  of  recent  advances  in  computational  aeroacousties  (CAA)  is  presented.  CAA  algorithms 
must  not  be  dispersive  and  dissipative.  It  should  propagate  waves  supported  by  the  Euler  equations  with 
the  correct  group  velocities.  Computation  domains  are  inevitably  finite  in  size.  To  avoid  the  reflection  of 
acoustic  and  other  outgoing  waves  at  the  boundaries  of  the  computation  domain,  it  is  required  that  special 
boundary  conditions  be  imposed  at  the  boundary  region.  These  boundary  conditions  either  absorb  all  the 
outgoing  waves  without  reflection  or  allow  the  waves  to  exit  smoothly.  High-order  schemes,  invariably, 
supports  spurious  short  waves.  These  spurious  waves  tend  to  pollute  the  numerical  solution.  They  must  be 
selectively  damped  or  filtered  out.  All  these  issues  and  relevant  computation  methods  are  briefly 
reviewed.  Jet  screech  tones  are  known  to  have  caused  structural  fatigue  in  military  combat  aircrafts. 
Numerical  simulation  of  the  jet  screech  phenomenon  is  presented  as  an  example  of  a successful 
application  of  CAA. 


1.  Introduction 

Aeroacoustic  problems  are  by  nature  very  different  from  standard  aerodynamics  and  fluid  mechanics 
problems.  Before  discussing  how  to  solve  aeroacousties  problems  numerically  or  simulate  them 
computationally,  an  approach  generally  referred  to  as  Computational  Aeroacousties  (CAA),  it  is 
important  to  recognize  and  to  have  a good  understanding  of  these  differences.  These  differences  pose  a 
number  of  major  challenges  to  CAA.  A few  of  the  important  computational  challenges  are  listed  below. 

a.  Aeroacousties  problems,  by  definition,  are  time  dependent,  whereas  aerodynamics  and  fluid  mechanics 
problems  are  generally  time  independent  or  involve  only  low  frequency  unsteadiness. 

b.  Aeroacousties  problems  typically  involve  frequency  range  that  spreads  over  a wide  bandwidth. 
Numerical  resolution  of  the  high  frequency  waves  becomes  a formidable  obstacle  to  accurate 

. numerical  simulation.  - 

c.  Acoustic  waves  usually  have  small  amplitudes.  They  are  very  small  compared  to  the  mean  flow. 
Oftentimes,  the  sound  intensity  is  five  to  six  orders  smaller.  To  compute  sound  waves  accurately,  a 
numerical  scheme  must  have  extremely  low  numerical  noise. 

d.  In  most  aeroacousties  problems,  interest  is  in  the  sound  waves  radiating  to  the  far  field.  This  requires  a 
solution  that  is  uniformly  valid  from  the  source  region  all  the  way  to  the  measurement-point  at  many 
acoustic  wavelengths  away.  Because  of  the  long  propagation  distance,  computational  aeroacousties 
schemes  must  have  minimal  numerical  dispersion  and  dissipation.  Also,  it  should  propagate  the  waves 
at  the  correct  wave  speeds  and  is  isotropic  irrespective  of  the  orientation  of  the  computation  mesh. 

e.  In  general,  flow  disturbances  in  aerodynamics  or  fluid  mechanics  problems  tend  to  decay 
exponentially  fast  away  from  a body  or  their  source  of  generation.  Acoustic  waves,  on  the  other  hand, 
decay  very  slowly  and  actually  reach  the  boundaries  of  a finite  computation  domain.  To  avoid  the 
reflection  of  outgoing  sound  waves  back  into  the  computation  domain,  radiation  boundary  conditions 
must  be  imposed  at  the  artificial  exterior  boundaries  to  assist  the  waves  to  exit  smoothly.  For  standard 
computation  fluid  dynamics  (CFD)  problems,  such  boundary  conditions  are  usually  not  required. 
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f,  Aeroacoustics  problems  are  archtypical  examples  of  multiple-scales  problems.  The  length  scale  of  the 
acoustic  source  is  usually  very  different  from  the  acoustic  wave  length.  That  is,  the  length  scale  of  the 
source  region  and  that  of  the  acoustic  far  field  region  can  be  vastly  different.  Computational 
aeroacoustic  methods  must  be  designed  to  deal  with  problems  with  greatly  different  length  scales  in 
different  parts  of  the  computational  domain, 

It  must  be  acknowledged  that  CFD  has  been  very  successful  in  solving  fluid  and  aerodynamics 
problems.  CFD  methods  are  generally  designed  for  computing  time  independent  solutions.  Because  of  the 
tremendous  success  of  CFD,  it  is  tempting  to  use  these  methods  to  solve  aeroacoustics  problems  as  well. 
In  the  past,  there  have  been  a number  of  attempts  to  do  just  that.  However,  the  results  have  proven  to  be 
quite  discouraging.  For  example,  Hsi  and  Perie  tried  to  use  a commercial  CFD  code  “RADIOSS”  to  solve 
sound  scattering  problems.  Their  effort  was  published  in  the  proceedings  of  the  Second  CAA  Workshop 
on  Benchmark  Problems  (NASA  CP3352,  June  1997,  edited  by  C.K.W.  Tam  and  J.C.  Hardin).  The 
results  were  disastrous.  The  computed  results  were  highly  dispersive  and  differed  significantly  from  the 
exact  solutions. 

As  elaborated  above,  it  should  be  clear  that  the  nature  of  aeroacoustics  problems  are  substantially 
different  from  those  of  traditional  fluid  dynamics  and  aerodynamics  problems.  To  be  able  to  compute  or 
simulate  aeroacoustics  problems  accurately  and  efficiently,  standard  CFD  schemes,  designed  for 
applications  to  fluid  problems,  are  generally  not  adequate.  For  this  reason,  there  is  a need  for  an 
independent  development  of  CAA.  This  was  what  happened  in  the  last  ten  years.  During  this  period,  great 
advances  have  been  made  both  in  the  development  of  CAA  methodology  and  in  their  applications  to  real 
world  problems. 

To  simulate  an  aeroacoustic  phenomenon  or  problem  numerically,  the  computation  algorithm 
must  consist  of  three  basic  elements.  They  are: 

(i)  A time  marching  computation  scheme. 

(ii)  An  artificial  selective  damping  algorithm  or  filtering  procedure. 

(iii)  A set  of  radiation/outflow  numerical  treatments  for  use  at  the  boundaries  of  the  computation  domain. 

A good  quality  time  marching  scheme  is  basic  to  any  computation  effort.  Artificial  selective 
damping  or  a filtering  procedure  is  essential  to  eliminating  spurious  numerical  waves  that  could 
contaminate  the  computed  solution.  Also  such  damping  terms  can  often  help  to  suppress  numerical 
instabilities  at  the  boundaries  of  the  computation  domain  or  at  surfaces  of  discontinuities  such  as  mesh- 
size-change  interfaces.  Numerical  boundary  treatments  serve  two  basic  purposes.  First,  they  are  to  allow 
outgoing  waves  to  leave  the  computation  domain  with  little  reflection.  Second,  they  are  to  reproduce  all 
the  effects  of  the  outside  world  on  the  computation  domain.  For  instance,  if  there  is  incoming  acoustic 
and  vorticity  wave  or  there  is  an  inflow,  they  are  to  be  generated  by  the  numerical  boundary  conditions. 

2.  Spatial  Discretization  in  Wavenumber  Space 


In  this  section,  we  consider  how  to  form  finite  difference  approximations  to  the  partial  derivatives  of 
spatial  coordinates.  The  standard  approach  of  truncated  Taylor  series  assumes  that  the  mesh  size  goes  to 
zero;  i.e.,  Ax->0,  in  formulating  finite  difference  approximation.  In  real  applications,  Av  is  never  equal  to 
zero.  Therefore,  we  will  investigate  how  well  we  can  approximate  3 Idx  by  a finite  difference  quotient 
with  a finite  Ax.  The  treatment  below  follows  the  work  of  Tam  and  Webb1  and  Tam2 
Suppose  we  use  a (2N+\)  point  stencil  to  approximate  3/3x;  i.e., 

^ ' -f„r  (2.1) 


(k)  i 


1 A' 

=t-  2>/m 


A Xy=_,V 


At  our  disposal  are  the  coefficients  a/,  j~N  to  N.  We  will  now  choose  these  coefficients  such  that  the 
Fourier  Transform  of  the  right  side  of  (2.1)  is  a good  approximation  of  that  of  the  left  side,  for  arbitrary 
function  f Fourier  Transform  is  defined  only  for  functions  of  a continuous  variable.  Let’s  first  generalize 
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(2.1)  to  all  points  along  the  x-axis.  The  generalized  form  of  (2.1),  applicable  to  any  set  of  points  at  Ax 
apart,  is 

1 N 

(*)  =~rHajf(x +Ax).  (2.2) 

AX 


(2. 1)  is  a special  case  of  (2.2).  By  setting  x = l Ax  in  (2.2),  (2.1)  is  recovered. 
The  Fourier  Transform  of  a function  F(x)  and  its  inverse  are  related  by 


F( a)  = — J F(x)e~,mdx , F(x)  = \~F(a)e,a'da . 


(2.3) 


Now  by  applying  Fourier  Transform  to  (2.2)  and  making  use  of  the  Derivative  and  Shifting 
Theorems,  it  is  found 


iaf  = 


Ax 


5>/“  l s=m 


j=-N 


where 


*=-r'Lvm' 

Ax 


(2.4) 


(2.5) 


Clearly,  "a  is  the  effective  wave  number  of  the  finite  difference  scheme,  a Ax  is  a periodic  function  of 
aAx  with  a period  of  2tc.  To  assure  that  the  Fourier  transform  of  the  finite  difference  scheme  is  a good 
approximation  of  that  of  the  partial  derivative  over  the  band  of  wavenumber  r|<aA  x<  tj,  it  is  required  that 
af  be  chosen  to  minimize  the  integrated  error,  E,  over  this  wavenumber  range,  where, 


E = j]«Ax -«Axf  d{aAx). 


For  a to  be  real,  the  coefficients  a,  must  be  antisymmetric;  i.e., 

«o  = 0, 

The  conditions  for  £ to  be  a minimum  are, 

dE  n 

da  ~ ’ 


a . — —a 

-i  j 


7 = 1,2, — ,N. 


(2.6) 


(2.7) 


(2.8) 


(2.8)  provides  N equations  for  the  A coefficients  a/,y=l,2,...,A'. 

For  a 7-point  stencil  (N= 3)  there  are  three  coefficients.  It  is  possible  to  require  (2.2)  to  be  accurate 
to  0(Ax4)  by  Taylor  series  expansion.  This  leaves  one  free  parameter  for  optimization.  By  setting  rpl.l, 
it  is  easy  to  find, 

a0  = 0.0,  t?,  = -GL,  = 0.77088238051822552, 

a2  = -a_2  = -0.166705904414580469,  a3  = -a,  = 0.02084314277031 176. 

The  choice  of  r\=l . 1 turns  out  to  provide  an  overall  best  group  velocity  approximation2. 

The  relationship  between  a Ax  and  otAx  over  the  interval  0 to  n for  the  optimized  7-point  stencil 
is  shown  graphically  in  Figure  1.  For  ccAx  up  to  about  1 .2,  the  curve  is  nearly  the  same  as  the  straight  line 
a Ax^oAx.  Thus,  the  finite  difference  scheme  can  provide  reasonable  approximation  to  ocAx<1.2  or 
A>5.2Ax  where  X is  the  wavelength.  For  aAx  greater  than  2,  the  a (a)  curve  deviates  increasingly  from 
the  straight-line  relationship.  Thus,  the  wave  propagation  characteristics  of  the  short  waves  of  the  finite 
difference  scheme  are  very  different  from  those  of  the  partial  differential  equations.  Figure  1 also  shows 
the  a Ax  versus  aAx  relations  for  the  standard  6th-order,  4,h-order  and  2nd-order  central  difference 
schemes.  Their  resolved  bandwidth  is  narrower.  In  other  words,  more  mesh  points  per  wavelength  is 
needed  for  propagating  a wave  accurately. 
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0.0  0.4  0.8  1.2  16  2.0  2.4  2.8  3.2 

aAx 

Figure  1.  a Ax  versus  a Ax.  optimized  ^-order  scheme, 

— * — • — ■ standard  6^-order  scheme,  — — standard  4^-crcirr  scheme.. 

standard  2n<*- order  scheme. 


Example.  Let  us  consider  the  initial  value  problem  of  the  one-dimensional  convective  wave  equation.  We 
will  use  dimensionless  variables  with  Ax  as  the  length  scale,  Axle  (c=  sound  speed)  as  the  time  scale.  The 
mathematical  problem  is: 

du  du 

-—■+  — = 0,  -oa<X<°°  - (2.9) 

at  ck 

/ = 0,  . u = f(x ) (2.10) 

The  exact  solution  is 

u = f(x-t). 

In  other  words,  the  solution  consists  of  the  initial  disturbance  propagating  to  the  right  at  a dimensionless 
speed  of  1 , without  any  distortion  of  the  waveform.  Let  us  consider  the  initial  disturbance  in  the  form  of  a 
Gaussian  function, 

frN2" 

/(jt)  = 0.5  exp 


(2.11) 


On  using  a (2N+l)-point  stencil  finite  difference  approximation  to  d/dx  of  (2.9),  the  original 
initial  value  problem  is  converted  to  the  following  system  of  ordinary  differential  equations 

du( 


dt 


+ YJaiua1=0’  (=  integer 


j=-N 

t = 0, 

For  the  Gaussian  pulse,  the  initial  condition  is, 


t = 0, 


tip  = 0.5  exp 


-in2\ 


(t2 

3 


(2.12) 

(2.13) 

(2.14) 
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The  system  of  ordinary  differential  equations  (2.12)  and  initial  condition  (2.14)  can  be  integrated 
in  time  numerically  using  the  Runge-Kutta  or  a multi-step  method.  The  solution  using  the  standard  2nd- 
order  (AM),  4th-order  (N- 2),  6th-order  (.M3)  and  the  optimized  7-point  scheme  (N=3f  T|=l.l)  at  f=400  are 
shown  in  Figure  2.  The  exact  solution  in  the  form  of  a Gaussian  pulse  is  shown  as  the  dotted  curves.  The 
second-order  solution  is  in  the  form  of  a wave  train.  There  is  no  resemblance  to  the  exact  solution.  The 
numerical  solution  is  totally  dispersed.  The  4th-order  solution  is  better.  The  6,h-order  scheme  is  even 
better.  But  there  are  still  some  dispersed  waves  trailing  behind  the  main  pulse.  The  optimized  scheme 
gives  a very  acceptable  numerical  result. 


central  difference  scheme,  (d)  DRP  scheme  (7-ooint  stencill. 
numerical  solution; exact  solution. 
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3.  Time  Discretization 


There  are  two  types  of  explicit  time-marching  schemes.  They  are: 


(1)  Single-step  scheme;  e.g.,  Runge-Kutta  method. 

(2)  Multi-step  scheme;  e.g.,  Adams-Bashford  method. 

Both  types  of  methods  are  discussed  in  standard  textbooks.  Here  we  will  only  discuss  the  optimized 
multi-step  method.  One  important  advantage  of  multi-step  method  is  that  it  can  be  used  in  a multi-size- 
mesh  multi-time-step  algorithm2.  Such  an  algorithm  is  very  efficient  and  nearly  optimal  in  computation 
time. 


Suppose  u(t)  is  the  unknown  vector.  The  time  axis  is  divided  into  a uniform  grid  with  time  step 
At.  We  will  assume  that  the  values  of  u and  du/dt  are  known  at  time  level  m,m-1,m-2,w-3.  (Note:  In  CAA, 
- du/dt  is  giveifby  the  governing  equations  of  motion.)  To  advance  to  the  next  time  level,  Tam  and  Webb1 
use  the  following  4-level  finite  difference  approximation 


u‘"+i)  = u'*,+A/£aJ 


.V »-» 


(3.1) 


The  last  term  on  the  right  side  of  (3.1)  may  be  regarded  as  a weighted  average  of  the  time  derivatives  at 
the  last  4 mesh  points.  There  are  four  constants,  namely,  h0,  hu  b2  and  />3.  They  can  be  determined  by 
optimization  through  the  use  of  Laplace  transform.  Tam  and  Webb  suggested  the  following  numerical 
values. 


ba  = 2.3025580888,  l\  =-2.4910075998, 
b2  = 1.5743409332,  b,  = -0.3858914222. 


4.  Radiation  and  Outflow  Boundary  Conditions 


Many  interesting  acoustic  problems  are  exterior  problems.  To  simulate  this  class  of  problems  it  is 
necessary  to  impose  radiation  and  outflow  boundary  conditions  at  the  boundaries  of  the  finite 
computation  domain.  To  ensure  that  the  computed  solutions  are  of  high  quality  these  boundary  conditions 
must  be  sufficiently  transparent  to  the  outgoing  disturbances  so  that  they  exit  the  computation  domain 
without  significant  reflections.  As  is  well  known,  the  linearized  Euler  equations  can  support  three  types  of 
waves.  Thus,  in  general,  the  outgoing  disturbances  would  contain  a combination  of  acoustic,  entropy  and 
vorticity  waves  each  having  a distinct  wave  propagation  characteristic.  Here  a set  of  radiation  and  outflow 
boundary  conditions  compatible  with  the  optimized  high-order  scheme  is  developed  starting  with  the 
asymptotic  solution.  A review  on  numerical  boundary  condition  can  be  found  in  Ref.  [4], 

Consider  the  exterior  problem  involving  a uniform  flow  of  velocity  u0  in  the  x-direction  and 
sound  speed  an  past  some  arbitrary  acoustic,  vorticity  and  entropy  sources  as  shown  in  Figure  3.  It  will  be 
assumed  that  the  boundaries  of  the  computation  domain  are  quite  far  from  the  sources.  At  boundaries 
where  there  are  only  outgoing  acoustic  waves  the  solution  is  given  by  the  following  asymptotic  solution, 

r i i 


(4.1) 
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where  F(0)=ao[M  cos  0 +(1  -M2  sin20)l/2].  The  subscript  ‘a’  in  (p a,ua,vapa)  above  indicates  that  the 
disturbances  are  associated  with  the  acoustic  waves  alone.  By  taking  the  time  t and  r derivatives  of  (4.1)  it 
is  straightforward  to  find  that  for  arbitrary  function  F of  the  acoustic  disturbances  satisfy  the  following 
equations. 


r i d { d _r 

, V(6)  dt  + dr  + 2r, 


u 

v 


( -S' 


= 0 + 0^ 


V 


(4.2) 


P 


Equation  (4.2)  provides  a set  of  radiation  boundary  conditions. 


outgoing 

vorticity 

and 

entropy 

waves 


Figure  3.  Acoustic,  vorticity  and  entropy  wave  sources  in  a uniform  flow. 


At  the  outflow  region,  the  outgoing  disturbances,  in  general,  consist  of  a combination  of  acoustic,  entropy 
and  vorticity  waves.  The  asymptotic  solutions  for  the  density,  velocity  and  pressure  fluctuations  are  given 

by 


X(x-u0t,y)  + pa 

p 

dw  . 

u 

~(x-uQt,y)  + ua 
dy 

V 

p 

l 

i 

L Pa  J 
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where  the  explicit  form  of  (p0,Ma,va,p„)  may  be  found  in  (4.1).  The  functions  y,  \j/  and  F are  entirely 
arbitrary.  These  functions  can  be  eliminated  by  combinations  of  derivatives.  In  this  way,  the  following 
outflow  boundary  conditions  can  be  derived. 


dp 

— + u, 
dt 


du 

dt 


■ + Un 


dp  __1 
° dx  a't 
du 
dx 
ch 


dt  ° dx 


dv  . 
dt  °dx 


1 dp 
Po  dx 
1 dp 


(4.4) 


Po  dy 


+ cos  + sin  + JL  - o 
V(0)  dt  dx  dy  2r 


5.  Artificial  Selective  Damping 


Numerical  waves  of  wavenumber  a for  which  a( a)  is  not  nearly  equal  to  a (see  figure  1)  will  not 
propagate  with  the  correct  wave  speed.  For  waves  on  the  right  side  of  the  maximum  of  the  a versus  a 
curve  of  figure  1,  the  wave  speed  is  negative  (because  da  Ida  is  negative)  or  opposite  to  the  correct  wave 
propagation  direction.  These  are  the  spurious  waves  of  the  computation  scheme5.  They  must  be  removed 
from  the  computation  if  a high  quality  numerical  solution  is  desired.  A way  to  do  this  is  to  apply  artificial 
selective  damping  to  the  discretized  equations.  Suppose  we  add  a damping  term,  D(x),  to  the  right  side  of 
the  momentum  equation  of  the  linearized  Euler  equations  in  one  dimension 


du 

dt 


^ 1 dp 
Po  dx 


= D(x). 


(5.1) 


Let  us  discretize  the  spatial  derivative  using  the  7-point  optimized  stencil. 

dut  1 'A  ^ V1  j (r 

<5'2) 

In  (5.2)  it  is  assumed  that  D(  is  proportional  to  the  values  of  wf  within  the  stencil  and  dj  be  the  weight 
coefficients  and  v is  the  artificial  kinematic  viscosity.  v/(Ac)2  has  the  dimension  of  (time)'1,  so  that  d/s 
are  pure  numbers.  Now  we  wish  to  choose  d/s  so  that  the  artificial  damping  would  be  concentrated 
mainly  on  the  high  wavenumber  or  short  waves. 

The  Fourier  transform  of  the  generalized  continuous  form  of  (5.2)  is, 

du  v ~ . _ 

-T+— =-77- (5.3) 
dt  (Ax) 

where 


D(cAx)  = ^ dj djl 


ijoAx 


(5.4) 


j=-  3 


On  ignoring  the  terms  not  shown  in  (5.3),  the  solution  is 


rD(aAx)l 

u~ew  . (5.5) 

Since  D{aAx)  depends  on  the  wavenumber,  the  damping  will  vary  with  wavenumber.  For  our  need,  we 
like  D to  be  zero  or  small  for  small  aAv  but  large  for  large  aAr.  This  can  be  done  by  choosing  dj 
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appropriately.  Tam,  Webb  and  Dong5  suggested  that  D should  be  a positive  even  function  of  a Ax.  They 
used  a Gaussian  function  as  a template  to  determine  d/s.  Their  analysis  gives  the  following  values. 

d{)  - 0.3276986608 

d]  = d_,  = —0.2357 1 88 1 5 

d2  =d_2  = 0.0861506696  ' 

di  =d^  =-0.0142811847 

A plot  of  D versus  ocAr  for  this  choice  is  shown  in  Figure  4.  Figure  4 indicates  that  the  grid-to-grid 
oscillations  with  ocAv=ji  or  wavelength  equal  to  2 mesh  spacings  is  most  heavily  damped.  On  the  other 
hand,  aAv=0  and  the  low  wavenumber  waves  are  hardly  damped  at  all.  Experience  has  shown  that  the 
adoption  of  artificial  selective  damping  in  a computation  scheme  can  effectively  eliminate  all  spurious 
short  waves  in  the  numerical  solution. 


Figure  4.  D(aAx)  versus  aAx  (damping  curve),  a = 0.37T,  7-point  stencil. 


3.2 
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6.  Numerical  Simulation  of  the  Jet  Screech  Phenomenon 

It  had  been  reported^7 * 9 * *  that  structural  fatigue  damages  sustained  by  the  Bl-A,  B2-B  and  the  F15  aircrafts 
were  most  likely  the  result  of  dynamic  loading  contributed  by  screech  tones.  The  screech  phenomenon  is 
known  to  occur  when  a supersonic  jet  is  imperfectly  expanded.  It  is  driven  by  a feedback  loop.  The  tones 
are  discrete  frequency  sound.  They  are  generated  by  the  interaction  of  downstream  propagating  instability 
waves  of  the  jet  flow  and  the  shock  cells  inside  the  jet  plume. 

Experimentally,  screech  tones  from  choked  jets  are  found  to  belong  to  four  distinct  modes.  The 
Ai  and  A2  modes  are  axisymmetric  modes  with  slightly  different  frequencies.  The  B and  C modes  are 
flapping  modes.  Their  frequencies  also  differ  only  slightly.  In  the  literature,  there  are  semi-empirical 
formulas  designed  to  predict  the  screech  frequencies.  However,  these  formulas  are  unable  to  predict 
individually  the  frequencies  of  the  Aj  and  A2  modes  nor  the  frequencies  of  the  B and  C modes.  They  give 
only  a single  approximate  frequency  for  the  axisymmetric  modes  and  a single  frequency  for  the  flapping 
modes.  The  feedback  loop  of  the  screech  tones  is  highly  nonlinear.  As  a result,  there  is  no  tone  intensity 
prediction  formula  available  at  the  present  time;  not  even  a totally  empirical  one. 

To  provide  numerical  prediction  of  the  screech  frequencies  and  intensities  as  well  as  to  better 
understand  the  physics  of  the  problem,  Shen  and  Tam10-12  performed  numerical  simulations  of  the  screech 
phenomenon.  For  marching  the  solution  in  time,  they  used  the  multi-size-mesh  multi-time-step 
Dispersion-Relation-Preserving  scheme.  This  is  an  improved  and  extended  version  of  the  optimized 
scheme  discussed  in  Sections  2 and  3.  Artificial  selective  damping  was  included  both  for  the  elimination 
of  spurious  short  waves  as  well  as  for  shock  capturing  purpose.  In  their  work,  a nonlinearized  form  of  the 
outflow  boundary  conditions  was  used.  Flow  entrainment  was  also  included  in  their  radiation  boundary 
condition.  To  simulate  the  effect  of  fine  scale  turbulence,  the  k-e  turbulence  model  was  employed. 
Screech  is  driven  by  a self-excited  feedback  loop.  Hence  no  external  excitation  is  needed.  In  the  computer 
simulation  of  Shen  and  Tam,  the  computation  started  from  zero  initial  condition.  In  other  words,  the  jet 
was  turned  on  at  time  zero.  After  the  transient  disturbances  had  exited  the  computation  domain,  the 
computation  locked  itself  onto  a screech  mode  (sometime  on  two  coexisting  modes)  automatically, 
exactly  as  in  a physical  experiment. 

Shen  and  Tamlft  12  reported  that  their  numerical  simulation  not  only  reproduced  the  screech 
frequencies  of  the  four  modes  correctly  but  also  the  tone  intensities  were  accurately  given  by  the 
computed  solution.  Experimentally,  it  is  observed  that  as  the  Mach  number  of  the  jet  increases,  at  certain 
jet  Mach  numbers,  there  is  an  abrupt  switch  of  the  screech  mode.  This  is  referred  to  as  the  staging 
phenomenon.  In  the  numerical  simulations  of  Shen  and  Tam,  the  staging  phenomenon  was  reproduced.  In 
addition,  they  studied  the  effect  of  jet  temperature  and  nozzle  lip  thickness  on  tone  frequencies  and 
intensities.  Their  computed  results  matched  well  with  those  measured  experimentally13.  Figures  5 and  6 
show  comparisons  between  the  screech  tone  wavelengths  and  intensities  measured  from  the  numerical 
simulations  of  Shen  and  Tam12  and  the  physical  experiments  of  Ponton  and  Seiner1’.  As  can  be  seen,  the 
CAA  simulations-  areJn  excellent  agreement  with  experimental  measurements  over  the  entire  range  of 
Mach  number  investigated. 


7.  Conclusion 

In  this  paper,  computational  challenges  for  the  solution  of  aeroacoustics  problems  are  discussed.  To  meet 
these  challenges,  specially  designed  high-order  CAA  computation  schemes  are  developed.  In  addition  to 
time  marching  schemes,  impressive  advances  have  also  been  made  on  numerical  boundary  treatment  as 

well  as  on  how  to  eliminate  spurious  numerical  waves.  As  an  example  to  illustrate  the  usefulness  and 

effectiveness  of  these  newly  developed  CAA  methods,  results  of  the  numerical  simulation  of  screech 
tones  from  under-expanded  supersonic  jets  are  briefly  discussed.  It  appears  that  CAA  has  now  emerged 

from  just  a novel  computation  method  to  become  a powerful  practical  tool.  CAA  will  undoubtedly  have 

an  enormous  impact  on  the  solution  and  investigation  of  practical  aeroacoustics  problems  in  the  years  to 

come. 
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Reference  # of  Paper:  Invited  Lecture  #1 
Discusser’s  Name:  Dr.  Brian  J.  Tester 

Author’s  Name:  Prof.  C.  K.  W.  Tam 

Question: 

The  numerical  simulations  show  energy  over  a broad  range  of  frequencies.  Is  the 
broadband  energy  meaningful  and,  if  so,  what  noise  generation  mechanisms  does  it  represent? 

Answer: 

Because  the  computation  took  a long  time  to  carry  out  the  simulation  was  not  as  long  as  I 
would  like.  The  relatively  short  simulation  time  causes  the  half-widths  of  the  screech  tones  in  the 
computed  noise  spectra  to  be  broader  than  they  should  be.  For  this  reason,  one  should  not  take 
the  broadband  part  of  the  spectra  too  seriously. 


Discusser’s  Name:  Prof.  Dr.  Rolf  Radespiel 

Author’s  Name:  Prof.  C.  K.  W.  Tam 

Question: 

By  using  a turbulence  model  in  the  acoustic  simulation  one  changes  the  nature  of  the 
unsteady  problem  to  be  solved.  Haw  can  you  be  sure  not  to  lose  relevant  physical  information  by 
introducing  a turbulence  model? 

Answer: 

In  the  numerical  simulation,  a k-e  turbulence  model  was  used  to  simulate  the  effect  of  the 
fine-scale  turbulence.  This  is  justified  because  of  scale  separation.  The  wavelengths  of  the 
instability  waves  in  the  jet  flow  are  much  larger  than  the  scale  of  the  fine-scale  turbulence.  This 
length  scale  disparity  permits  one  to  ignore  the  individual  eddy  motions  of  the  fine-scale 
turbulence.  Only  the  collective  effect  of  the  fine-scale  turbulence  on  the  large-scale  motion  is 
considered. 


Discusser’s  Name:  Dr.  Mahmood  Khalid 

Author’s  Name:  Prof.  C.  K.  W.  Tam 

Question: 

Did  you  have  any  difficulties  with  computational  instabilities  at  the  interface  between  the 
Euler  and  Navier-Stokes  domains?  Did  you  apply  the  selective  artificial  damping  technique 
uniformly  across  the  entire  computational  field? 

Answer: 

No,  I have  never  encountered  computational  instabilities  at  the  interface  between  Euler 
and  Navier-Stokes  regions.  In  the  numerical  simulation  of  jet  screech,  background  artificial 
selective  damping  was  incorporated  in  the  computation  in  each  computational  subdomain. 


